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Abstract (Basic) : EP 624847 A 

The generator produces an approximation function based on pairs of 
values (Xi,Yi) associating a dependent magnitude (Yi) with an 
independent magnitude (Xi) , and generates new data pairs (Xa,Ya) 
following the approximation. 

The generator has a first part (10) that iteratively determines a 
current linear regression function to select those of the current 
regression functions that deliver the approximation to the current set 
of input values with least error. The parameters (p,q) of the best 
regression function are then retained, and used by the second part (17) 
to generate the new approximate pair of data values. 

ADVANTAGE - Requires little computational capacity and gives rapid 
results under controlled maximum error, for use in neural networks or 
in approximation to measured values in industrial processes. 
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Abstract (Equivalent) : US 5519647 A 

An apparatus (5) for generating an approximation function based on 
first pairs ({XI, Yl) to (X6, Y6) ) of values associating a dependent 
variable (Yl to Y6) with an independent variable (XI to X6) , and for 
determining second pairs (XA, Y'A) of values of said variables in 
accordance with said approximation function, characterized in that the 
apparatus comprises : 

first means (10) : 

for iteratively determining at least one current linear regression 
function by making first errors of alternating sign equal in absolute 
value, which first errors have been measured between, respectively, 
first values (Y3, Y4 , Y5) of the dependent variable for three pairs 
(X3, Y3) (X4, Y4) (X5, Y5) of a series of pairs, and second values 
(Y'3, Y'4, Y'5) of the dependent variable determined, in accordance 
with said current linear regression function, for the same values (X3, 
X4, X5) of the independent variable, 

for selecting that one of the current linear regression functions 
which produces an approximation of all the pairs of said series with a 
minimum error, and 

for coding the selected linear regression function with the aid of 
specific codes (p, q) , and 

second means (17) for determining said second pairs (XA, Y'A) with 



the aid of said specific codes. 
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© Dispositif et methode pour generer une fonction d'approximation. 

© Dispositif (5) pour generer une fonction d'approximation fondee sur des premiers couples ((X, , Y t ) a (X*. 
Y 6 )) de valeurs associant une grandeur dependante (Yi a Y 6 ) a une grandeur independante (X, a Xg), et pour 
determiner des seconds couples (X A , Y' A ) de valeurs desdites grandeurs d'apres ladite fonction d'approximation. 
Le dispositif comprend : 

- des premiers moyens (10) pour determiner iterativemant au moins une fonction lineaire courante de 
regression, pour selectionner celle des fonctions lineaires courantes qui delivre ('approximation de tous les 
couples de ladite suite avec des erreurs minimales. et pour coder, a I'aide de codes specifiques (p, q), la 
fonction lineaire de regression selectionnee, 

- et des seconds moyens (17) pour determiner lesdits seconds couples (X Al Y' A ) a I'aide desdits codes 
specifiques. 

Le dispositif peut etre utilise pour le calcul de valeurs approchees de fonctions mathematiques, par exemple 
dans un reseau de neurones, ou pour determiner une fonction de regression formant une approximation de 
valeurs experimentales de mesure, par exemple des mesures dispersees issues du suivi d'un processus 
industriel. 

L'invention concerne egalement une methode pour generer une fonction d'approximation. 
Application : calcul approche de fonctions mathematiques, calcul de fonctions de regression. 
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L'mvention concerne un dispositif et une methode pour generer une fonction d'approximation. ce!te-ci 
etant fondee sur des premiers couples de valeurs associant une grandeur dependante a une grandeur 
independante, et pour determiner des seconds couples de valeurs desdites grandeurs a partir de ladite 
fonction d'approximation. 

5 Un dispositif et une methode de ce type sont connus du brevet US-A- 3 789 203 qui decrit un 

generateur de fonctions operant une approximation par interpolation iterative. Ce dispositif est prevu pour 
etre utilise dans des applications de traitement de donnees necessitant un caicul de fonctions telles que 
sin(x), tg(x) par exempie. Ce dispositif ne requiert qu'une capacite de stockage minimale de la part d'un 
dispositif utilisateur. A partir de deux points appartenant a une fonction a interpoler, la methode tout d'abord 
to interpote la fonction par une droite reliant les deux points, puis fait une approximation des ecarts entre la 
droite et la fonction par des approximations polynomials de degre croissant. Ensuite elle substitue aux 
points initiaux des points approximatifs pour reduire la longueur du segment reliant les points a trailer et 
enfin reitere les operations precedentes. 

Une telle methode necessite des ressources importantes en moyens de caicul et ne peut etre mise en 
75 oeuvre qu'avec des calculateurs performants. 

Or, il existe des applications ou une teMe methode ne peut pas etre mise en oeuvre car elles ne 
disposent pas des ressources suffisantes. De plus, pour certaines applications on peut se satisfaire d'un 
caicul approche de la fonction pour des valeurs, en nombre limite, de la grandeur independante. 

II peut s'agir d'une fonction sigmoVde appliquee a des potentiels neuronaux delivres par au moins un 
20 neurone dans un reseau de neurones. II peut s'agir d'une autre fonction non lineaire, par exempie une 
fonction racine, pour des calculs de distances entre des etats de neurones. Les applications peuvent aussi 
concerner d'autres dispositifs comme des generateurs de fonctions, des calculateurs ou autres. 

Pour calculer une telle fonction, sans passer par une fonction d'approximation, on peut utiliser 
differentes manieres. 

25 On peut effectuer le calcui mathematique exact pour chaque valeur de la grandeur independante a 
traiter, en programmant un calculateur selon les methodes connues. Une telle methode necessite d'effec- 
tuer a chaque fois les memes operations ce qui peut necessiter beaucoup de temps si le nombre de 
valeurs est eleve. 

On peut aussi prealablement stocker dans une memoire des tables precalculees. Dans ce cas, la 
30 lecture en memoire du resultat peut §tre rapide. Mais pour couvrir avec un pas assez fin toutes les valeurs 
possibles de la grandeur independante, il faut alors disposer de tables de grandes capacites. Ces methodes 
de caicul presentent done des inconvenients. 

D'autre part, on peut etre conduit a identifier deux grandeurs qui sont dans la dependance Tune de 
Tautre par des couples de valeurs associant une grandeur dependante a une grandeur independante. Ainsi. 
35 dans le suivi d'un processus industriel, on peut etre conduit a mesurer par exempie un rendement R d'une 
operation en fonction de la temperature T a laqueile a ete realisee ladite opeation R = f(T). Un suivi du 
processus peut consister a consigner des lots de couples de mesures sur un graphique. Celui-ci peut etre 
utilise pour caracteriser le processus ou pour en deduire de nouveaux parametres de reglage de ladite 
operation. Ceci fait par exempie I'objet de I'article de H. ISHIBUCHI et H. TANAKA, "Regression analysis 
40 with interval model by neural networks" in "IEEE International Joint Conference on Neural Networks", vol. 2, 
18-21 Nov. 1991, SINGAPORE. II est souhaitable que ces nouveaux parametres soient representatifs du 
fondement de ladite operation et que les fluctuations de mesures inherentes a ce genre de processus 
soient ecartees. It est done souhaitable de determiner une approximation de la fonction f(.)- 

Ainsi dans un cas il peut s'agir de mesures erratiques ou entachees d'erreurs que Ton desire 
45 representer par une fonction d'approximation. 

Dans un autre cas, on connait des valeurs precises mais I'utilisation a en faire ne necessite pas une 
grande precision et une fonction d'approximation suffit. 

Un des buts de ('invention est de generer une fonction d'approximation avec des moyens materiels 
reduits permettant de calculer rapidement un nombre reduit de valeurs de la grandeur dependante utiles a 
so I'application sans avoir pour cela a determiner d'autres valeurs de la fonction d'approximation. Un but 
complementaire est de delivrer des valeurs qui peuvent etre approchees dans la limite d'une erreur 
maximale controlee. 

Ce but est atteint avec un dispositif caracterise en ce qu'il comprend : 
- des premiers moyens : 

55 - pour determiner iterativemant au moins une fonction lineaire courante de regression en rendant 

egales, en valeur absolue, des premieres erreurs de signes alternes mesurees entre des premieres 
valeurs de la grandeur dependante pour trois couples d'une suite desdits premiers couples, et 
respectivement des secondes valeurs de ta grandeur dependante determinees, d'apres ladite 
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fonction lineaire courante. pour ies mimes vaieurs de la grandeur independante, 

- pour selectionner celle des fonctions lineaires courantes qui delivre ('approximation de tous Ies 
couples de ladite suite avec des erreurs minimales. 

- et pour coder, a I'aide de codes specifiques, la fonction lineaire de regression selectionnee, 

- et des seconds moyens pour determiner lesdits seconds couples a I'aide desdits codes specifiques. 
Ainsi avantageusement on determine une fonction lineaire de regression approchant au mieux Ies 

differents couples de vaieurs connus. Les resultats approches ainsi delivres torment un compromis 
satisfaisant pour de nombreuses utilisations du dispositif generateur de fonctio... 

Une fonction lineaire de regression est une fonction simplificatrice qui represente un phenomene 
complexe en reduisant les parametres significatifs. En representant la suite de couples de vaieurs par des 
points dans un espace a deux dimensions, la fonction lineaire de regression devient une droite de 
regression. 

Ainsi apres avoir defini la droite de regression par des codes, on peut calculer une valeur approchee de 
la grandeur dependante en tout point de la droite de regression avec des moyens reduits pour des vaieurs 
quelconques de la grandeur independante. 

L'invention concerne egalement une methode pour generer une fonction d'approximation, la methode 
comprenant : 

- une premiere phase : 

- pour determiner iterativement au moins une fonction lineaire courante de regression en rendant 
egales. en valeur absolue, des premieres erreurs, de signes alternes, mesurees entre des 
premieres vaieurs de la grandeur dependante pour trois couples d'une suite desdits premiers 
couples, et respectivement des secondes vaieurs de la grandeur dependante determinees. d'apres 
ladite fonction lineaire, pour les memes vaieurs de la grandeur independante, 

- t pour selectionner celle des fonctions lineaires courantes qui delivre I'approximation de tous les 
couples de ladite suite avec des erreurs minimales, 

- et pour coder la fonction lineaire de regression selectionnee a I'aide de codes specifiques, 

- et une seconde phase pour determiner lesdits seconds couples a I'aide desdits codes specifiques. 
Les moyens mis en oeuvre par l'invention peuvent etre formes par un calculateur programme ou par un 

circuit dedie. lis peuvent aussi mettre en oeuvre des neurones. 

Un dispositif mettant en oeuvre des neurones selon l'invention peut etre utilise par un reseau de 
neurones, dont il peut notamment en constituer un sous-ensemble. En effet, pour fonctionner, le reseau de 
neurones doit disposer de moyens pour appliquer une fonction non lineaire d'activation aux potentiels de 
neurones qu'il delivre. Selon l'invention, le dispositif muni de neurones peut calculer une approximation de 
cette fonction non-lineaire d'activation. II peut egalement calculer des distances entre des etats de neurones 
en calculant une approximation d'une fonction racine carree destinee a etre exploitee dans le reseau de 
neurones. 

Lorsque la taille de la suite de couples de vaieurs fournis initialement est elevee, on peut diviser la suite 
de couples en plusieurs sous-ensembles pour determiner plusieurs droites de regression et ameliorer la 
precision de I'approximation. La fonction d'approximation de la suite de couples est alors formee par une 
fonction lineaire par morceaux pour laquelle une exigeance de continuite entre les morceaux peut etre ou 
non imposee. 

Certains couples de I'ensemble de couples de vaieurs peuvent avoir une influence particuliere que Ton 
peut concretiser en donnant un coefficient de ponderation specifique a chaque couple. Dans ce cas, I'erreur 
affectee a chaque couple tient compte de ce coefficient de ponderation specifique. 

Ces differents aspects de l'invention et d'autres encore seront apparents et elucides a partir des modes 
de realisation decrits ci-apres. 

L'invention sera mieux comprise a I'aide des figures suivantes donnees a titre d'exemples non limitatifs 
qui represented : 

Figure 1 : un graphique montrant une representation a deux dimensions d'un ensemble de points avec 
une droite de regression D. 

Figure 2 : un graphique montrant un ensemble de points et des droites servant a la determination d'une 
enveloppe. 

Figure 3 : un organigramme d'une premiere variante de mise en oeuvre de la methode a partir de 
triplets de points. 

Figure 4 : un organigramme d'une seconde variante de mise en oeuvre de la methode a partir de triplets 
de points. 

Figure 5 : une partie d'organigramme d'une troisieme variante de mise en oeuvre de la methode a partir 
de couples de points. 
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Figure 6 : une partie d' organigramme indiquant une preselection de points appartenant a une envefoppe 
inferieure ou a une enveloppe superieure de I'ensemble des points. 

Figure 7 : un organigramme pour la determination des enveloppes inferieure et superieure. 
Figure 8 : un schema d'un dispositif selon ('invention. 
5 F,gure 9 : un schema d'un dispositif de codage pour le calcul et le codage dune droite de regression 
dans le cas general. 

Figure 10 : un schema d'un dispositif de codage pour le calcul et le codage dune droite de regression 
pour la premiere variante. 

Figure 11 : un schema d'un dispositif de transcodage pour le calcul des valeurs de la grandeur 
w dependante a partir d'un ensemble de droites de regression codees. 

Figure 12 : un schema d f une realisation neuronale des moyens de com'paraison de la figure 10. 
Figure 13 : une representation montrant une fonction d'approximation formee de plusieurs droites. 
Figure 14 : deux representations montrant une determination des limites de validite de deux droites 
consecutives. 

75 Figure 15 : une representation concernant le raccordement de droites consecutives. 

L'invention concerne I'approximation d'une fonction connue uniquement a travers un certain nombre de 
points par exemple P, ... P 6 (figure 1) dans une representation a "deux dimensions. Chacun de ces points 
est defini par un couple de valeurs (x, y) reliant la grandeur independante x a la grandeur dependante y. 
Par la suite, il sera question d'une paire (respectivement d'un triplet) de points, ce qui met en oeuvre deux 
20 couples (respectivement trois couples) de valeurs. On ordonne les points d'apres un ordre croissant des 
valeurs d'abscisses X h ce qui definit un indice i croissant avec lesdites valeurs. Une convention inverse 
peut etre faite en modifiant en consequence i'expose ci-apres. 

Selon l'invention on effectue une approximation de I'ensemble des couples (Xi, Y,), (X 2 , Y 2 ) par une 
droite de regression D ayant pour equation : 
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(1) D \ y - p.x + q. 



ou x et y sont des variables courantes. 

Pour cela on considere dans une premiere variante trois couples de valeurs, par exemple (X 3) Y 3 ), (X*. 
Ya). (X 5 , Ys), et on determine une droite de regression D par equilibrage des ^reurs absolues. Une erreur 
est mesuree par la difference apparaissant, pour une abscisse x donnee, entre la valeur y du point et 
I'abscisse y mesuree sur la droite de regression. Equilibrer les erreurs sur trois points, consiste a avoir trois 
erreurs egales en valeur absolue avec un signe d'erreur oppose aux deux autres pour le point ayant une 
abscisse x comprise entre les abscisses x des deux autres points. Puis on examine si pour les points 
35 restants de I'ensemble, I'erreur qui les separe de la droite reste inferieure ou egale en valeur absolue a 
I'erreur prealablement determined pour les trois points selectionnes. Ceci decoule du fait que I'on 
s'interesse a une erreur p.re cas relative a I'ensemble de tous les points a prendre en consideration c'est- 
a-dire la plus grande erreur, en valeur absolue. qui existe entre un des points et la droite de regression Si 
toutes les erreurs sont effectivement inferieures ou egales. la droite est selectionnee pour representer les 
points s.non on recommence les operations avec trois nouveaux couples de valeurs pour determiner une 
autre droite de regression. 

^ II peut exister plusieurs droites de regression representant tous les points de I'ensemble Selon la 
methode on determine la droite de regression optimale qui minimise I'erreur pire cas definie prealablement 
La figure 1 represente un exemple forme de six points Pi a P 6 disposes selon une representation a 
deux dimensions. A des fins duplication considerons le resultat final. On observe que la droite de 
regression D de la figure 1 est situee de sorte que les erreurs sont egales en valeur absolue pour les points 
P 3 . P* et P 5 . Pour les points P, , P 2 , P 6 les erreurs sont inferieures aux precedentes. en valeur absolue 
Dans un ensemble de points P, a P 6 , la methode va ainsi consister a rechercher les trois points 
parlicuhers, ici P 3 , P., P 5 , qui permettent de determiner la droite de regression optimale qui minimise 
I erreur p ,re cas puis a coder cette droite. Dans le cas de la figure 1 representant un resultat final si on 
trace deux droites D, et D 2 , paralleles a la droite de regression D, qui passent respectivement par les 
points P 3 . P 5 , d'une part et P. d'autre part, on constate que tous les points de I'ensemble sont a I'interieur 
d un bandeau hmite par les droites D, et D 2 ou sur ces droites. 

La phase de determination de la droite de regression peut donner lieu a plusieurs mises en oeuvre dont 
55 seules les plus avantageuses seront decrites ci-apres. 
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Premiere variante de la premiere phase de la methode . 

La figure 3 represente ia suite des etapes a mettre en oeuvre pour determiner la droite de regression 
optimale. a 

5 Parm. ('ensemble des po.nts P,... P N (bloc 100). on selectronne (bloc 102) trois points quelconques P, 

P,-. P k avec i < j < k. Ces trois points servent a determiner la droite de regression qui minimise I'erreur sur 
y pour ces tro.s po.nts. Cette determination est effectuee de maniere analytique preferentiellement par des 
moyens programmes. On determine une droite de recession D pour que trois erreurs (bloc 104) • 
E PD (P ; , D), E P0 (P r D), Ep D (P kf D) 

;o entre la droite D et chacun des points venfient : 

E PD (P. D) = -Epo (Pj, D) = Ep D (P k ,D) 
avec 

Epo (Pi, D) = Y; - (p.Xj + q) 

et des relations analogues pour les autres erreurs. 

On determine la droite de regression D a I'aide des coefficients p et q de I'equation (1) tels que : 
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p = — - I 



et 
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= x i- Yi : X i- Y 3 + Y _i - x :-y k * x k . y, - x f . y k 

2(X k - Xj) " 



L'erreur se rapportant a un triplet (Pi, Pj, P k ) s'ecrit alors : 
35 E T (P„ P h P k ) = | Epo (P„ D) | 

Lorsque I'erreur E T a ete ainsi calculee, on examine si les autres points de I'ensemble engendrent des 

nlTp 'tr 6 "^ 0 " ^ a ' eS ' e " Val6Ur abS °' ue - 6 Ce " eS des points P " P i- p - Pour cela on selectionne un 
M I 1 h , C 106) * ddi,ionnel et on calcule ,a va| eur absolue de I'erreur E Pm (bloc 108) apparaissant entre 
40 la valeur de la variable y au point P m et la droite D. 

Lorsque cette erreur E Pm est inferieure ou egale en valeur absolue a E T (bloc 110) (repere Y) le point 
P m addit,onnel est accep.e et la methode se poursuit (bloc 112) avec un point additionnel suivant (bloc 
106 S, tous les po.nts additionnels satisfont au critere | E Pm | s E T . la droite D est accep.ee et ses 
coeffic.ents sont utilises pour coder la droite de regression optimale D opl (bloc 114) 
-5 Lorsque cette erreur E Pm est superieure a E T (bloc 110) (repere N). le triplet de points P, P, P k 

T™T e ,u, nof aCCeP -' 6t 3Utre ,riP ' et de P ° in,S (Dl0C 116) est Choisi < lien 101 > dans ''ensemble 
de points (bloc 102). La methode s'acheve a I'obtention de la droite D (bloc 114, sa.isfaisant ce critere 
meme si tous les triplets de points n'ont pas ete examines. 

50 e, J, l ,™ t rt aPP ! ra T? S !' Uati0nS 00 ^ '' iSSUe de ''^ ape 116 tous les P° in,s P° ssiDles °«* ete examines 
so et aucun des tnplets na fourni de solution (bloc 118). Dans ce cas. il est possible de reprendre le 
deroulement de la premiere variante en remplacant le test du bloc 1 10 par le test suivant : 

\E Pm \S a.E T 

ou a est un coefficient legerement superieur a 1. Dans cette meme situation, il est egalement possible de 
taire appel a la deuxieme variante. 

La premiere phase peut comprendre les etapes suivantes : 
A - selection de trois couples de valeurs parmi ladite suite. 
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B - caicul de la fonction lineaire courante de regression D et determination d'une erreur de triplet 

E T -\E PD \ s'y rapportant, 
C - selection d'un couple additionnel, 

D - calcul d'une erreur additionnelle E Pm entre le couple additionnel et ladite fonction. 
E - et lorsque | E Pm | < E^HO) pour le couple additionnel, la methode reprend a I'etape C avec un 

couple additionnel suivant, 
F - et lorsque | E Pm \ > E T pour au moins un couple additionnel, 

la methode reprend a I'etape A avec une selection d'un autre groupe de trois couples parmi ladite suite 
de couples, 

G - et lorsque | E Pm \ < E T pour tous les couples additionnels, la fonction lineaire courante de 

regression est codee et stockee en tant que fonction lineaire d'approximation. 
On peut choisir de scruter ('ensemble des triplets en prenant un ordre croissant ou un ordre decroissant 
ou un ordre aleatoire pour effectuer cet examen. Le triplet qui sera retenu pour determiner la droite de 
regression pourra de ce fait etre detecte a un instant quelconque du deroulement de cette scrutation. II 
75 s'ensuit que la rapidite d'obtention de la droite de regression depend de I'instant au cours duquel le triplet 
est detecte. Sa mise en oeuvre presente un degre de complexite allant de N a N 4 ou N est le nombre de 
points initiaux. Sa complexite est done reduite pour un petit nombre de points. Cette variante permet 
d'obtenir une realisation materielle avec une forte parallelisation. Elle est tres peu sensible a une troncature 
des valeurs et fournit un resultat exact. 
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Deuxieme variante de la premiere phase de la methode. 



Dans cette deuxieme variante (figure 4), on selectionne successivement des triplets de points, on 
calcule chaque fois une droite de regression, et, par recurrence, on selectionne celle qui delivre I'erreur E P0 
la plus grande e'est-a-dire correspondant a I'erreur pire cas pour I'ensemble de points considered. 
Selon la seconde variante, la premiere phase comprend les etapes suivantes : 
A - selection de trois couples de valeurs parmi ladite suite. 

B - calcul d'une fonction lineaire courante de regression D et determination d'une erreur de triplet 
£ T =|£p D | s'y rapportant. 

^ c " comparaison de I'erreur E T avec une erreur optimale E op initi^see a une valeur strictement 
negative, 

D - et lorsque E T > E op , mise a jour de I'erreur optimale E op en remplacant E op par E T et mise a jour 
des codes d'une fonction lineaire optimale de regression D op en remplacant ceux-ci par les 
codes de la fonction lineaire courante D, 
5 E - puis retour a I'etape A pour selectionner trois autres couples, 

F - et lorsque tous les triplets de couples de valeurs de la suite ont ete testes, les derniers codes de 

la fonction lineaire optimale D op constituent les codes de la fonction lineaire d'approximation. 
^ A chaque examen d'un triplet, on compare I'erreur E T dudit triplet avec I'erreur optimale precedemment 
memorised et on met a jour I'erreur optimale E op avec la plus grande valeur de I'erreur E T determined pour 
> chaque triplet. On met a jour egalement les parametres de la droite optimale. Avant la mise en oeuvre de 
I'etape A, il faut initialiser la valeur E op a une valeur faible negative, par exemple -1. 

Dans ce cas, la rapidite d'obtention de la droite de regression ne depend pas du mode de scrutation 
des triplets. Sa complexite est de degre N 3 . Une mise en oeuvre materielle peut beneficier de la grande 
regularity de I'algorithme mis en oeuvre. II est peu sensible aux troncatures des donnees et fournit une 
solution exacte. 

Troisieme variante de la premiere phase de la methode. 

Dans cette troisieme variante (figure 5), on selectionne d'abord une paire de points a laquelle on ajoute 
un point supplemental, situe entre ces deux points, afin de former un triplet de points. Pour cela. on 
modifie les etapes A, B et C de la premiere variante. les autres etapes restant les memes. Les etapes 
modifiees sont telles que : 

A1 - modifie I'etape A en operant une selection de deux couples de valeurs appartenant a ladite 
suite, tel qu'il existe au moins un couple additionnel intermediate ayant une grandeur indepen- 
dante (X, - Xe) comprise entre les grandeurs independantes dudit couple pour constituer au 
moins un triplet de couples, 

A2 - modifie I'etape A premierement en determinant une fonction lineaire annexe qui contient les 
deux couples selectionnes et deuxiemement en determinant des secondes erreurs entre les 
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granaeurs dependantes des couples intermediates possibles et ladite fonction imeaire annexe : 
et lorsque ces secondes erreurs sont toutes de meme signe. selection du couple 
intermediate fourntssant la plus grande seconde erreur, en valeur absoiue. pour former un 
triplet de couples de valeurs forme du couple intermediate et des deux couples 
selectionnes. 

. et lorsque ces secondes erreurs sont de signes differents. reprise de la methode a I'etape 
A1, 

B1 - I'etape B est effectuee avec ledit triplet selectionne, 

C1 - modifie I'etape C en selectionnant un couple additionnel dont la grandeur independante n'est 
pas comprise entre la grandeur independante des deux couples selectionnes. 

Lorsque I'erreur | E Pm | est superieure a I'erreur E T (bloc 110). la methode reprend a I'etape A1 (lien 
101) avec une nouvelle selection de paire de points (bloc 102a). 

On observe que le deroulement de cette troisieme variante depend de la scrutation des valeurs et done 
des valeurs elles-memes. La complexite de la mise en oeuvre de cette variante varie entre N et N 3 ce qui 
lui confere un certain avantage par rapport aux variantes precedentes. Une forte parallelisation des moyens 
de mise en oeuvre peut etre operee mais I'implementation des moyens peut presenter un manque de 
regularity ce qui peut constituer un handicap pour une realisation integree. Cette variante est peu sensible 
aux erreurs d'arrondis des valeurs et fournit une solution exacte. 

Quatrieme variante de la premiere phase de la methode. 

Elle concerne la determination de la droite de regression a partir des enveloppes. 

II est possible de reduire le nombre de triplets a examiner en determinant des enveloppes respective- 
ment superieure et inferieure entourant les points extremes dans la representation bidimensionnelle de 
I'ensemble de points. Une enveloppe superieure ou une enveloppe inferieure est definie telle qu'en joignant 
par une droite deux points adjacents quelconques de I'enveloppe, tous les autres points soient situes d'un 
meme cote respectivement de I'enveloppe superieure ou de I'enveloppe inferieure. On determine ainsi tous 
les points appartenant a ces dites enveloppes. 

La determination de la droite de regression va consister a considerer les paires de points adjacents 
d'une des enveloppes auxquels on associe un point intermediate n'appartenant pas a ladite enveloppe 
pour constituer un triplet et operer comme cela vient d'etre decrit dans le cas des paires de points de la 
troisieme variante. Si une solution optimale n'a pas ete trouvee, on considere les paires de points de i'autre 
enveloppe. 

Pour mettre en oeuvre une enveloppe, on selectionne une paire de points adjacents appartenant a 
I'enveloppe. On determine alors s'il existe un point dispose de telle fagon que son abscisse soit 
intermediate entre les abscisses des points selectionnes. Lorsque ce point n'existe pas on passe a une 
autre paire de points de la meme enveloppe. Pour certaines paires, lorsqu'il apparait qu'il existe un ou 
plusieurs de ces points intermediates, on choisit le point intermediate le plus eloigne de la droite 
contenant les deux points de la paire pour former un triplet et pour determiner une droite de regression. 
Pour determiner si cette droite de regression peut etre selectionnee comme droite de regression optimale. 
la methode met en oeuvre les memes operations que celles decrites prealablement dans le cas de la 
troisieme variante. 

Pour cela on modifie la troisieme variante telle que (figure 6). prealablement a I'etape 102a (figure 5), la 
premiere phase de la methode comprend une etape (bloc 100a) de determination d'une enveloppe 
inferieure et'ou d'une enveloppe superieure reunissant les points les plus extremes de I'ensemble de 
points, la selection des paires de points a I'etape 102a etant faite parmi les points adjacents appartenant a 
Tune ou I'autre enveloppe. La selection de ladite paire de points est effectuee lorsqu'il existe au moins un 
point intermediate ayant une abscisse situee entre les abscisses des points de la paire de points. S'il 
existe plusieurs points intermediates, on forme le triplet avec le point intermediate le plus eloigne de la 
droite passant par les deux points qui forment la paire de points. Si une solution n'est pas trouvee avec la 
premiere enveloppe. on poursuit le traitement avec la seconde enveloppe. 

La determination des enveloppes inferieure et superieure est effectuee selon I'organigramme de fa 
figure 7. L'indice des points allant croissant avec la variable x, le premier point P 0 fait ainsi partie des deux 
enveloppes. Considerons d'abord I'enveloppe inferieure. les points appartenant a I'enveloppe inferieure sont 
reperes par la lettre Q. Un point courant Q v est repere par l'indice v. 

A Petape 400 les deux premiers points de I'enveloppe sont : Q c = P 0 et Q, = P, . Le dernier point 
courant Q v est repere : v = 1 pour P, . Le nombre k de points traites est comptabilise. 
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Etape 402 : un premier test est effectue pour determiner si le dernier point P N _, a ete traite pour 
detecter !a fin de ia determination de I'enveloppe. 

Etape 404 : dans ie cas contraire on teste si v > 1. Lorsque v< 1 on incremente v (v = v + 1) et on 
prend le point P k comme point Q v (etape 407). On incremente I'indice k pour trailer le point P suivant 
(etape 409) et reprise de la methode a I'etape 402. 

Lorsque v > 1, on calcule une droite passant par les points Q v _ , et Q v (etape 406) et on determine le 
signe c de I'erreur sur la variable dependante entre le point P k courant et cette droite (Q v _,, Q v ). Ceci a 
pour but de determiner si le point courant est au-dessus ou en dessous de la droite (Q v . Q v ). 

Lorsque le signe t £ o. il faut supprimer le dernier point Q v . decrementer v tel que v = v - 1 (etape 
410) et reprendre la methode a I'etape 404. On peut de cette maniere etre conduit a supprimer certains 
points deja acceptes lorsqu'un point suivant necessite de les supprimer. 

Lorsque le signe « est strictement positif, la methode reprend a I'etape 407 avec un point suivant. 

Cet organigramme s'applique a la determination des enveloppes inferieure et superieure en inversant te 
signe de I'erreur a considerer. 

Pour faire comprendre les mecanismes ainsi mis en oeuvre considerons, a titre d'exemple, le cas 
simple forme par les points Pi, P 2 , P 3 . Pa de la figure 2 et determinons I'enveloppe inferieure. Le point P: 
est le premier point de I'enveloppe d'ou Q 0 = Pi . On prend P 2 comme deuxieme point Q1 . On calcule la 
droite L 0 reliant Pi, P 2 . Le segment Pi, P 2 est accepte comme premier segment de I'enveloppe basse. On 
examine si P 3 est au-dessus de la droite L 0 , ce qui est le cas sur la figure 2. Le segment P 2 P 3 est accepte 
comme deuxieme segment de I'enveloppe basse. On calcule ensuite la droite Li reliant P 2 et P 3 et on 
examine si P 4 est au-dessus de la droite Li . Ce n'est pas le cas sur la figure 2. Ceci entraTne que le point 
P 3 et la droite Li sont d'abord ecartes. De meme, P 4 est en dessous de la droite L G done le point P 2 et la 
droite L D sont ensuite aussi ecartes comme ne pouvant pas appartenir a I'enveloppe inferieure et les 
operations se poursuivent avec la determination de la droite L 2 reliant P,, P«, puis I'examen des points 
suivants (non represented). 

La complexite de mise en oeuvre de cette variante basee sur des enveloppes est de degre N 2 et 
depend de I'ordre de scrutation des donnees. Cette complexite est moindre que celle des variantes 
precedentes et de ce fait fournit un resultat rapidement. La regularity de ^implementation est moyenne mats 
cette variante est peu sensible aux erreurs d'arrondi des valeurs et delivre une solution exacte. 

Pour certaines applications. »' peut etre souhaitable d'accroitre la precision de la determination de la 
fonction d'approximation dans certains domaines de la grandeur independante x et d'affecter des coeffi- 
cients de ponderation W { aux points P,. C'est par exemple le cas lorsque la fonction d'approximation est 
faiblement variable avec la grandeur independante x. A certains points peuvent alors etre affectes des 
coefficients de ponderation. Ceux-ci peuvent etre communs a plusieurs points ou etre individuels pour 
chaque point. Ces coefficients de ponderation Wj sont, par la suite, considered strictement positifs. 

Dans les cas ou des coefficients de ponderation existent, on definit une erreur E PD entre un point Pj et 
la droite de regression D tel que : 

E PD (P h D) = W,.[Y t -{pX + q)] 

ou E PD est une valeur signee. La determination de la droite de regression D pour trois points Pj, P jt P k est 
alors modifiee en ce que les quantttes p et q de I'equation (l).deviennent : 

(2) p = NUMP/DET et q = NUMQ/DET 

ou les quantites NUMP, NUMQ et DET sont definies par : 

DET = W i .W i .(X r X i ) + Wj.WK.fXk-Xj) + W t .\N k .(X k -X,) 

NUMP = Wi.W^YpYi) + Wj.W k .(Y k -Yi) + W t .W k .(Y k -Y t ) 

NUMQ = Wi.Wj.tXj.YrXj.Yj) + ^.W^y^y^ + Wi.W k .(X k .Y r X,Y k ) 

Par ailleurs. I'erreur E T associee a ce triplet peut etre exprimee et calculee par : 
avec 

EjiPi,P h P h ) = \E PD (P it D\ 

Epd(P„D) = -E PD (Pj,D) = E PD (P k ,D) 

ou 
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DET 



La premiere et la seconde variante de la methode de la premiere phase decrites precedemment 
peuvent etre mises en oeuvre en faisant intervenir les coefficients de ponderation ci-dessus. Cette mise en 
oeuvr^ peut etre effectuee en programmant un calculateur. 

Dans le cas ou l*on selectionne une paire de points en faisant intervenir des coefficients de ponderation 
W, la methode presente les amenagements suivants. 

Pour les points Pj. P,, P k qui forment le triplet, on definit un facteur F, qui lie les coordonnees 
respectives et les coefficients de ponderation respectifs des points P i( P,. P k . Au point central 1, on associe 
le facteur F, tel que : 




Ce facteur F, intervient dans la determination du point intermediate P, qui sera a selectionner (i < I < k) 
pour constituer le triplet en tenant -compte des poids affectant les points. On choisit un point P, et on 
calcule la droite de regression Dt associee au triplet P„ P, t P k et I'erreur E T i associee au triplet. 

Cette variante de la methode consiste a modifier uniquement I'etape A2 de la troisieme variante (bloc 
102b, figure 5). Cette etape determine Texistence et la valeur d'un point intermediaire servant a former un 
triplet. On cherche d'abord a former une droite de regression situee en dessous de P, et de P k . Pour 
chaque point intermediaire P,, on determine si F, = 1 et E P0 (Pj, D,) < 0. Si au moins un point verifie cette 
condition, il n'existe pas de droite de regression situee en dessous de Pj et P k> sinon on determine une 
grandeur G max qui est la quantite maximale parmi : 

. d'une part les quantites E P0 (P it D p ) pour tous les points intermediates, 

. d'autre part les quantites (a) suivantes, uniquement pour les points intermediaires pour lesquels F, < 1 
avec 



S'il existe au moins un point intermediaire tel que F, > 1. on definit egalement une grandeur G min qui 
est la valeur minimale prise par (a) uniquement pour les points intermediaires tels que F, > 1. On examine 
alors s'il existe au moins un point intermediaire tel que : 

Ep D (P jt D ( ) > 0 

G max ~ Or E T | 

et G min . a > E T » (a, coefficient £ 1). S'il existe un tel point, il est choisi comme point intermediaire pour 
former le triplet. 

S'il n'existe pas de point intermediaire tel que F t > 1, on teste s'il existe au moins un point' 
intermediaire tel que : 
Epd (Pi, 0\) > 0 

G max < a. E T r (a : coefficient > 1). 

Si un tel point existe, il est choisi comme point intermediaire pour former le triplet. 

Si aucun triplet n'a ete constitue, on cherche a former une droite de regression situee au-dessus des 
points Pj et P k . La meme methode est reprise en inversant le signe des erreurs E PD . 

Si aucun point P, ne peut etre selectionne. on recommence avec une autre paire P it P k . 

Lorsque I'ensembie de points a traiter est trop important pour etre represents par une seule droite de 
regression, la methode code alors plusieurs droites de regression chacune etant determined selon la 
methode decrite precedemment. 

La figure 13 represente un exemple dans lequel la fonction d'approximation est constitute de plusieurs 
droites de regression. 
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Dans une premiere situation, de par la connaissance des grandeurs a trailer, on peut vouioir imposer 
des limites a chaque droite. de regression suivant la grandeur independante x. Ainsi on peut vouioir 
disposer d'une droite de regression D 0 entre les valeurs [X a ,X b [ de la grandeur, borne x a incluse. borne x b 
exclue. De meme avec D b pour [X b , X c [ et D c pour [X c , X^ . Dans ce cas le probieme revient a determiner 
5 une droite dans un domaine limite et a appliquer a chaque fois la methode deja decrite. 

Mais il est possible de faire que pour chaque droite, la methode determine des limites suivant la 
grandeur independante x sans qu'elles soient fixees au depart. 

Le principe de la determination d'une limite optimale entre deux droits de regression adjacentes est 
represents sur la figure 14-A. Soient deux droites de regression D1 et D2 non optimales. La droite D1 est 
w determined a partir de N1 points et la droite D2 est determinee a partir des N2 points restants avec N1 + 
N2 = N, ou N est le nombre total de points. Les droites D1 et D2 forment une approximation avec 
respectivement des erreurs maximales E1 et E2. Supposons que E2 soit superieure a E1 et que Ton 
cherche a reduire E2, il faut alors transferer vers D1 un point appartenant a D2. Lorsque le nombre de 
points diminue, I'erreur resultante est soit inchangee soit reduite. Lorsque le nombre de points augmente, 
is I'erreur resultante est soit inchangee soit augmentee. En fonction du nombre n de points appartenant aux 
droites D1 et D2, les variations des erreurs E1 et E2 peuvent presenter les allures representees sur la figure 
14-B pour I'exemple choisi. II s'ensuit que I'approximation obtenue globalement par D1 et D2 est optimale 
lorsque ces erreurs sont voisines, c'est-a-dire situees dans la zone delimitee par un trait discontinu sur la 
figure 14-B. Ces erreurs ne sont pas necessairement egales compte tenu du caractere discontinu des 
20 valeurs de la grandeur independante. 

Pour determiner la valeur d'abscisse X Um entre les deux droites : 

- on determine la droite D1 sur un certain nombre de points et on calcule I'erreur E1 maximale, 

- on determine la droite D2 sur les points restants et on calcule I'erreur E2 maximale, 

- on compare E1 et E2 et on transfere un point de la droite qui presente la plus forte erreur vers la 
25 droite qui presente la plus faible erreur, 

- on determine la valeur limite X, im lorsqu'il se produit une inversion dans le rapport entre les deux 
erreurs. 

Cette methode s'applique par recurrence a plusieurs droites delimiters entre elles par plusieurs limites 
d'abscisses. 

30 Les droites sont determinees a partir d'une suite discrete et limitee de grandeurs de mesure. 
Neanmoins, pour I'exploitation des droites de regression, il est necessaire de definir leur domaine 
d'existence qui s'etend sur un continum de valeurs situees entre deux limites d'abscisses. Or la determina- 
tion des droites fournit une suite de droites qui ne se raccordent pas necessairement par leurs extremites. II 
peut apparaitre utile pour certaines applications d'eviter qu'un ecart apparaisse sur la valeur de la grandeur 

35 dependante y pour des valeurs voisines (X Mm - f ) et (X| im + *) de la grandeur independante ou € est une 
tres faible valeur. II est possible de choisir que I'abscisse limite X, im appartienne exclusivement a Tune ou a 
I'autre droite. II est aussi possible de faire suivre la determination des suites de droite d'une procedure de 
raccordement de droites. 

Ceci est represents sur fa figure 15. Une solution peut consister a valider la droite D1 jusqu'a I'abscisse 
40 du premier point appartenant a la droite D2 et a recalculer la droite D2 a partir de la nouvelle valeur de 
grandeur dependante y ainsi determinee. En conservant le dernier point appartenant a D2, on determine 
ainsi une nouvelle droite D'2 representee en pointille sur la figure 15. Une procedure analogue peut etre 
appliquee pour substituer la droite D'3 a la droite D3. On obtient ainsi un ensemble forme de plusieurs 
droites de regression forme, dans I'exemple decrit, par les droites D1, D'2, D'3. Cet ensemble constitue une 
45 approximation des grandeurs de mesure en reduisant une erreur maximale entre les grandeurs de mesure 
et ('ensemble de droites. Une variante plus adaptee a cette reduction consiste a faire que ce soient les 
droites correspondant aux plus fortes erreurs qui imposent leurs points limites comme points limites aux 
droites correspondant a des erreurs plus faibles. 

50 Dispositif cour une mise en oeuvre de la premiere phase 

Un schema d'un dispositif pour generer une fonction d'approximation selon I'invention est represents 
sur la figure 8. Le dispositif 5 regoit des couples (Xj, Y t ) associant la grandeur dependante Yj a la grandeur 
independante X 4 . Les couples (X it Yj) entrent dans les premiers moyens 10 pour determiner et coder la 
55 fonction lineaire de regression formant une approximation des couples de mesure. Les codes specifiques 
ainsi determines sont transmis (connexion 9) aux seconds moyens 17 qui determinent des seconds couples 
(X A , Y* A ) d'apres les codes specifiques a la fonction d'approximation. 
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On distingue d'une part la seconde vanante qui determine toutes les droites de regression passant par 
chaque connbinaison realisable de triplets avec la suite de premiers couples, et les autres variantes pour 
iesquelles a chaque droite de regression determmee <assoaee a son erreur E T ), on examine si les points 
additionnels restants fournissent bien des erreurs addttionnelles inferieures ou egales a I'erreur de triplet E T . 

La figure 9 represente un dispositif adapte a la seconde variante. II comprend : 

- une memoire 12c MEM qui stocke notamment tous les points appartenant a {'ensemble de points a 
traiter. Ces points sont represents par leurs coordonnees (x, y), et eventuellement leurs poids W ou 
ieurs inverses 1/W, 

- une unite de calcul 13c COMPUT qui calcule, pour chaque triplet selectionne, la droite de regression 
adaptee a chaque triplet, c'est-a-dire les codes p. q de la droite et I'erreur E T associee au triplet. 

De plus, un controleur 11c CONTR permet d'effectuer la gestion des operations et d'adresser des 
nouveaux triplets en effectuant les lectures/ecritures de la memoire 12c et le chargement de l'unite de 
calcul 13c par de nouveaux triplets. La selection de la droite de regression qui est a conserver, c'est-a-dire 
pour cette variante celle delivrant la plus grande erreur de triplet, est effectuee par I'unite de calcul 13c. 

La figure 10 correspond au cas des autres variantes pour Iesquelles pour chaque droite courante de 
regression on examine si les autres points additionnels de la suite delivrent une erreur inferieure a I'erreur 
de triplet. 

Pour cela, les premiers moyens 10 comprennent une unite de calcul 13c et une unite de comparaison 
14c formant des moyens de calcul 19c. 

L'unite de calcul 13c. transmet les codes p, q de la droite de regression du triplet courant a I'unite de 
comparaison 14c COMPAR qui determine si les points additionnels de I'ensemble de points generent une 
erreur plus faible avec cette droite de regression que celle generee par les points du triplet courant. Pour 
cela, l'unite de comparaison 14c effectue le test : 

I E Pm \> E T 

Si la droite de regression courante n'est pas acceptee (test positif), un autre triplet est selectionne et une 
procedure analogue est a nouveau effectuee. Si le test est negatif pour tous les points additionnels, la 
droite de regression est acceptee et ses parametres sont charges par I'unite de calcul 13c dans la memoire 
12c. 

Pour les variantes qui operent a partir d'une paire de points, l'unite de calcul 13c est programmee pour 
determiner les points intermediates et pour former les triplets issus des paires de points, ceux-ci etant eux- 
memes eventuellement issus d'enveloppes de points. Dans ce dernier cas, I'unite de calcul 13c est 
egalement programmee pour determiner les enveloppes. L'unite de comparaison 14c effectue ensuite les 
comparaisons des erreurs se rapportant aux differents triplets. 

Pour effectuer la gestion des operations et adresser des nouveaux triplets (lorsque des triplets 
examines n'ont pas donne satisfaction), le dispositif comprend, comme precedemment, un controleur 11c 
qui gere : 

- les lectures/ecritures de la memoire 12c, 

- le chargement de l'unite de calcul 13c par de nouveaux triplets, 

- I'appel vers l'unite de comparaison 14c de tous les points additionnels a examiner successivement. 

II est possible de mettre en oeuvre la premiere phase de la methode a I'aide de moyens de 
comparaison 14c (figure 12) ayant une structure neuronale. Considerons le cas le plus complet ou il existe 
des coefficients de ponderation associes a chaque point P t . Lorsque ces coefficients n'existent pas il suffit 
de !eur donner une valeur unite dans les explications qui vont suivre. A titre d'exemple, on peut stocker 
dans la memoire 12c (figures 9, 10) les inverses des coefficients Wj affectes individuellement aux points P,. 
L'unite de calcul 13c determine selon les equations (2) et (3), pour un triplet donne, les codes : -p, -q, E T et 
■ E T . 

La condition a tester dans I'unite 14c pour un point additionnel P m a tester, de parametres X m , Y m , W m , 

est : 

Y m - pX m - (Er'W m ) - q > 0 
ou 



Y m - p.X m + (E T /W m ) - q < 0. 
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Ces tests sont aisement realisabies par une unite 14c ayant une organisation neuronale. En effet, ce 
test est constitue de fonctions lineaires. de fonctions de seuiliage et de fonctions logiques aisement 
realisees dans un reseau de neurones. 

Une telle unite neuronale de comparaison est representee sur ta figure 12. Elle comprend trois 
neurones N1, N2, N3. Les neurones N1 et N2 recoivent des donnees E1, E2. E3. E4. Le neurone N3 regoit 
les sorties des neurones N1 et N2. Chacune des entrees de ces neurones est affectee d'un coefficient 
synaptique Ci selon la technique connue mise en oeuvre dans les reseaux de neurones. Cette technique 
est par exemple decrite dans R.P. LIF?MANN, "An introduction to computing with neural nets" IEEE ASSP 
Magazine, April 1987, pp. 4 to 22. 

Pour effectuer les tests precedents, on programme les neurones N1 et N2 selon le tableau I : 

TABLEAU I 



Neurone 


E, 


Ci 


E 2 


c 2 


E 3 


c 3 


E* 


C* 


1 


Y m 


1 


X m 


- p 


1 w m 


-E T 


1 


- Q 


2 


Y m 


1 


X m 


- p 


1A/V m 


E t 


1 


- Q 



Tous les coefficients synaptiques du neurone N3 sont egaux a 1. 
Les neurones N1 et N2 calculent chacun un potentiel neuronal s tel que : s = E, Q.E, 
A ce potentiel s, il faut appliquer une fonction d'activation A. Pour les neurones N1 et N3, la fonction 
d'activation Ai est : 

At (S) = \ Si S > 0 
A,(s) = 0 si s < 0. 



Pour le neurone N2, la fonction d'activation A 2 est : 



A 2 (s) = 1 si s < 0 



A 2 (s) ■= 0 si s k 0. 



La sortie du neurone N3 vaut 1 si la condition a tester est satisfaite et vaut 0 dans le cas inverse. 

On observe sur le tableau I que les donnees -p, -q, -E T .+ E T sont des codes qui apparaissent comme 
coefficients synaptiques dans le schema de la figure 12. Les donnees Y m , X m , 1/W m sont des donnees qui 
apparaissent sur les entrees du meme schema. Lorsque tous les points additionnels ont ete testes et 
qu'une dro.te de regression a ete selectionnee les codes p, q, W charges dans la memoir© 12c sont alors 
utilises pour mettre en oeuvre la seconde phase de la methode (figure 11). 

L'interet presente par la realisation neuronale decrite est de pouvoir paralleliser les differentes 
operations a mettre en oeuvre selon les variantes deja decrites. Le fonctionnement d'une telle realisation 
neuronale est alors tres rapide. 

Lorsque les droites de regression ont ete determines, la seconde phase consistant a calculer des 
seconds couples (X A , Y' A ) de valeurs des grandeurs est alors mise en oeuvre dans des moyens de 
decodage 17 (figure 11). Les codes des droites sont charges dans une memoire 12a qui au cours de la 
seconde phase est adressee. par un controleur 11 a, pour fournir les codes des droites de regression 
adressees. La memoire 12a, organisee en lignes par exemple, contient pour chaque droite de regression 
les parametres p, q et x L ou x L est la limite superieure d'abscisse pour laquelle chaque droite de reqression 
est definie. y 

La memoire 12a contient ainsi un tableau de parametres 



X L1' p o' q o 

x Lm' p m-r q m-l 

correspondant aux m droites de regression stockees. 



12 
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Les moyens oe decodage 17 comprennent : 

- le controleur 1 1 a , 

- la msmoire 12 0 , 

- et une unite de decodage I3 a . Celle-ci regoit une valeur de requete X A de la grandeur independante 
5 pour laqueile on desire obtenir un resultat Y' A issu de la fonction d'approximation. Pour cela le 

controleur 1 1 a presente successivement certaines iignes de codes de la memoire 12 a et, pour chaque 
ligne de codes. I'unite de decodage I3 a teste le code x L pour determiner si la valeur X A de la 
grandeur d'entree peut etre traitee par la droite de regression codee par ladite ligne. Si cela n'est pas 
le cas, I'unite 13 3 requiert du controleur 1l a une autre ligne de codes. Lorsque cette condition est 
;o verifiee. I'unite 13 a effectue le calcul : 

Y' A = p (V X A + q n 

si la ligne d'ordre n est la ligne adequate. 
75 Comme cela vient d'etre decrit, les premiers moyens 10 peuvent constituer un dispositif de codage et 
les seconds moyens 17 peuvent constituer un dispositif de decodage, les deux dispositi f s etant separes. lis 
peuvent en effet etre distants run de I'autre, les operations "de codage etant separees des operations de 
decodage. Mais les deux moyens 10 et 17 peuvent aussi former un dispositif de codage/decodage unique. 
Dans ce cas certains organes peuvent etre affectes successivement a plusieurs operations. Notamment : 
20 - il peut s'agir d'une memoire 12 unique qui peut stocker les parametres des couples de valeurs et les 
codes des droites de regression determinees. Les unites 12a et 12c forment alors la memoire 12, 

- il peut s'agir d'un controleur 1 1 unique qui regroupe les fonctions des contrdleurs 1 1c et 11a, 

- I'unite de calcul 13c et I'unite de decodage 13a peuvent egalement former une seule unite 13. 

II est egalement possible de realiser les moyens de decodage 17 (figure 11) selon une organisation a 

25 structure parallele. Une solution consiste alors a comparer, en parallele, la grandeur de requete X A a tous 
les codes x L . On obtient un ensemble de signaux T L representant les resultats de ces tests de comparaison. 
On realise ensuite, en parallele, pour toutes les valeurs de L, un test qui combine deux a deux les signaux 
T L et T L+1 et qui indique si la grandeur de requete X A est dans I'intervalle d'indice L Ce test est verifie pour 
une unique valeur de L, notee n. On accede alors aux codes p n et q n correspondants, et on en deduit la 

30 grandeur Y' A de la meme maniere que precedemment. 

Les premiers moyens 10 (de codage) associes aux seconds moyens 17 (de decodage) peuvent etre 
utilises pour determiner la valeur d'une fonction d'approximation par au moins une droite de regression. 
Ces determinations peuvent etre effectuees pour des valeurs quelconques de la grandeur independante (a 
I'interieur de bornes predefines definissant I'etendue de Taction de chaque droite de regression). Cette 

35 methode evite de stocker inutilement des tables de valeurs, pour lesquelles toutes les valeurs ne seront pas 
exploitees. Selon la methode on ne determine que les seules valeurs necessaires a I'application. La 
methode selon I'invention presente I'interet de ne calculer que les valeurs necessaires. Le dispositif de 
I'invention peut etre utilise en combinaison avec un processeur neuronal pour calculer une approximation 
d'une fonction non lineaire par exemple une fonction sigmoTde. Un processeur neuronal calcule des 

40 potentiels neuronaux issus de sommes de produits de coefficients synaptiques par des etats de neurones. 
La technique est developpee dans le document de R.P. LIPPMANN deja cite. Chaque potentiel neuronal 
doit etre soumis a I'action d'une fonction non lineaire. Pour cela, chaque potentiel neuronal constitue la 
grandeur independante decrite precedemment et le calcul de la fonction th(x) par exemple est effectue pour 
les seules valeurs de requetes utiles. La methode fournit des valeurs approximees, mais cette approxima- 

45 tion peut etre obtenue avec une precision qui peut etre accrue selon le nombre de droites de regression 
calculees et eventuellement les coefficients de ponderation W. La precision reste neanmoins tributaire de la 
precision des coordonnees des couples de depart. 

Une telle methode est particulierement interessante pour les calculs de fonctions connues (telles que 
des fonctions mathematiques) ou de fonctions inconnues explicitement, par exemple une fonction represen- 

50 tee par des points de mesures, que I'on desire simplifier par une fonction lineaire de regression. 

Dans sa version neuronale, I'invention est interessante dans des applications neuronales car elle 
apporte non seulement un traitement homogene mais aussi une grande compacite de ('architecture 
materielle necessaire. 

II est possible de combiner les variantes de la methode, dans leurs mises en oeuvre materielles pour 
55 combiner leurs avantages. Ainsi par exemple, les combinaisons suivantes sont avantageuses : 

Dans le cas ou les poids existent, deux combinaisons apparaissent interessantes : 
Dans la premiere combinaison, on utilise tout d'abord la premiere variante basee sur des triplets de 
points. Si une solution n'est pas trouvee. on enchaTne avec la deuxieme variante basee aussi sur des 
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triplets de pomts. Ceci permet d'obtenir les avantages suivants : 
Quelles que soient les donnees. on obtient une solution. 

Dans le cas ou on exige qu'on obtienne toujours une solution exacte. cette combinaison permet 
d'obtenir la rapidite maximale. En effet. dans le cas general ia solution est fournie directement par 
la premiere variante appelee, qui est celle fournissant le plus rapidement une solution exacte dans 
le cas ou les poids existent. La deuxieme variante est beaucoup plus lente mais elle permet de 
garantir que la combinaison fournit toujours une solution. 
Dans la deuxieme combinaiscn, on utilise tout d'abord la variante basee sur des paires de points. Si 
elle ne trouve pas une solution (qui peut etre approchee), on utilise la deuxieme variante basee sur 
des triplets de points. Ceci permet d'obtenir les avantages suivants : 

Quelles que soient les donnees. on obtient une solution (eventuellement approchee). 
Cette combinaison est en moyenne plus rapide que la precedente. 
En contrepartie, elle ne fournit generalement qu'une solution approchee. 
Dans le cas ou les poids n'existent pas, une combinaison interessante consiste a utiliser tout d'abord la 
variante basee sur I'enveloppe des points. Si elle ne trouve pas de solution, on utilise la deuxieme variante 
basee sur des triplets de points. Ceci permet d'obtenir les avantages suivants : 
Quelles que soient les donnees, on obtient une solution, " ■ 
Cette combinaison permet d'obtenir la meilleure rapidite moyenne. 
Lorsque la fonction lineaire de regression a ete determinee selon la methode decrite, il est possible de 
determiner non seulement des grandeurs Y' A ressortissant a ladite fonction mais egalement des grandeurs 
ressortissant a la derivee de ladite fonction. En effet. pour une valeur de requete X A de la grandeur 
independante, il suffit d'extraire des moyens de memorisation la valeur p ressortissant a la valeur X A pour 
connaftre la valeur de cette derivee correspondante. 

Revendications 

1. Dispositif (5) pour generer une fonction d'approximation fondee sur des premiers couples ((Xi , Yi) a 
(Xe, Y e )) de valeurs associant une grandeur dependante (Y, a Y 5 ) a une grandeur independante (Xi a 
Xs), et pour determiner des seconds couples (X At Y* A ) de valeurs desdites grandeurs d'apres ladite 
fonction d'approximation caracterise en ce que le dispositif comprend : 

- des premiers moyens (10) : 

- pour determiner iterativemant au moins une fonction lineaire courante de regression en 
rendant egales, en valeur absolue, des premieres erreurs, de signes alternes, mesurees entre 
des premieres valeurs (Y 3 , Y 4 , Y 5 ) de la grandeur dependante pour trois couples (X 3t Y 3 ) (X*. 
Y*) (Xs, Y 5 ) d'une suite desdits premiers couples, et respectivement des secondes valeurs 
(Y' 3 . Y'a, Y' 5 ) de la grandeur dependante determinees, d'apres ladite fonction lineaire courante. 
pour les memes valeurs (X 3 , X*. X 5 ) de la grandeur independante, 

- pour selectionner celle des fonctions lineaires courantes qui delivre ('approximation de tous les 
couples de ladite suite avec des erreurs minimales, 

- et pour coder, a I'aide de codes specifiques (p, q), la fonction lineaire de regression 
selectionnee, 

- et des seconds moyens (17) pour determiner lesdits seconds couples (X A , Y' A ) a I'aide desdits 
codes specifiques. 

2. Dispositif selon la revendication 1 caracterise en ce que les premiers moyens (10) comprennent : 

- des moyens (12c) pour memoriser les grandeurs des premiers couples et lesdits codes specifi- 
ques, 

- des moyens (13c) (19c) pour calculer successivement les codes specifiques (p, q) de fonctions 
lineaires courantes de regression et leurs premieres erreurs, egalisees en valeur absolue, 
respectives (E T ) et pour selectionner au moins une fonction lineaire courante de regression, 

- et des moyens (11c) pour controler les moyens (12c) pour memoriser et les moyens (13c) pour 
calculer. 

3. Dispositif (5) selon les revendications 1 ou 2 caracterise en ce que pour selectionner ladite fonction 
lineaire, les premiers moyens comprennent des moyens (14c) : 

- pour determiner pareillement des erreurs additionnelles pour les autres couples (Xi . Y,) (X 2 . Y 2 ) 
(Xg, Y & ) de valeurs de ladite suite, 

- pour comparer les erreurs additionnelles avec lesdites premieres erreurs egalisees (E T ), 
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- et pour selectionner la fonction lineaire courante qui fournit des erreurs additionnelles inferieures 
ou egales en vaieur absoiue auxdites premieres erreurs egalisees (E T ). 

Dispositif selon la revendication 2 caracterise en ce que. pour selectionner la fonction lineaire de 
regression, les moyens pour calculer (13c) selectionnent celie des fonctions lineaires courantes de 
regression qui delivre des premieres erreurs egalisees (E T ) maximales. 

Dispositif selon une des revendications 1 a 4 caracterise en ce que les seconds moyens (17) 
comprennent : 

- des moyens (12a) pour memoriser lesdits codes specifiques d'au moins une fonction lineaire de 
regression, 

- des moyens (13a) pour selectionner, sur requite d'une vaieur de grandeur independante (X A ), 
une fonction lineaire de regression ressortissant a ladite vaieur, et pour calculer, a partir des 
codes specifiques regus des moyens pour memoriser, une vaieur de grandeur dependante (Y' A ) 
associee a la vaieur de grandeur independante (X A ), 

- et des moyens (1 1 a) pour adresser aux moyens pour calculer, et cur leur demande, les codes 
specifiques des fonctions lineaires de regression. 

Dans un dispositif pour generer une fonction d'approximation fondee sur des premiers couples ((Xj , Yi ) 
- (Xg, Y 6 )) de valeurs associant une grandeur dependante {Y^ - Y 6 ) a une grandeur independante (Xi - 
Xs), methode pour generer ladite fonction d'approximation a partir des premiers couples (X ; , Y t ) (X&, 
Y 6 ) de valeurs, et pour determiner des seconds couples (X Af Y' A ) de valeurs desdites grandeurs 
d'apres ladite fonction d'approximation, caracterisee en ce que, la methode comprend : 
,- une premiere phase : 

- pour determiner iterativement au moins une fonction lineaire courante de regression en 
rendant egales, en vaieur absoiue, des premieres erreurs E P0 , de signes alternes, mesurees 
entre des premieres valeurs (Y 3 , Y*. Y 5 ) de la grandeur dependante pour trois couples (X 3 , 
Y3), (X4, Y*) (X 5 , Y 5 ) d'une suite desdits premiers couples, et respectivement des secondes 
valeurs (Y' 3> Y* 4| Y's) de la grandeur dependante determined, d'apres ladite fonction lineaire, 
pour les memes valeurs (X 3 , X*, X 5 ) de la grandeur independante, 

- pour selectionner celle des fonctions lineaires courantes qui delivre I'approximation de tous les 
couples de ladite suite avec des erreurs minimales, 

- et pour coder la fonction lineaire de regression selectionnee a I'aide de codes specifiques, 

- et une seconde phase pour determiner lesdits seconds couples (X A , Y' A ) a I'aide desdits codes 
specifiques. 

Methode selon la revendication 6 caracterisee en ce que la premiere phase comprend les etapes 
suivantes : 

A - selection (102) de trois couples (P3, 94. P5) de valeurs parmi ladite suite, 

B - calcul (104) de la fonction lineaire courante de regression D et determination d'une erreur de 

triplet £Y = |£pd| s'y rapportant, 
C - selection (106) d'un couple additionnel (P1), (P2), (P6), 

D - calcul (108) d'une erreur additionnelle E Pm entre le couple additionnel et ladite fonction, 

E - et lorsque | E Pm | < £t<110) pour le couple additionnel, la methode reprend a I'etape C avec 

un couple additionnel suivant. 
F - et lorsque | E Pm \ > E T pour au moins un couple additionnel (1 10), la methode reprend a 

I'etape A avec une selection d'un autre groupe de trois couples parmi ladite suite de 

couples, 

G - et lorsque | E Pm | < E T pour tous les couples additionnels (112), la fonction lineaire courante 
de regression est codee et stockee (1 14) en tant que fonction lineaire d'approximation. 

Methode selon la revendication 6 caracterisee en ce que la premiere phase comprend les etapes 
suivantes : 

A - selection (102) de trois couples (P3, P4, P5) de valeurs parmi ladite suite. 

B - calcul (104) d'une fonction lineaire courante de regression D et determination d'une erreur de 

triplet E T = \E PC \ s'y rapportant, 
C- comparaison (210) de I'erreur E T avec une erreur optimale E op initialised a une vaieur 

strictement negative, 
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D - et iorsque E T > E op . mise a jour (l 14) de i'erreur optimale E op en remplagant E op par E T et 

mise a jour des codes d'une fonction lineaire optimale de regression D op en remplagant 

ceux-ci par les codes de la fonction lineaire courante (D), 
E - puis retour (1 16) a I'etape A pour selectionner trois autres couples (102), 
5 F - et Iorsque tous les triplets de couples de valeurs de la suite ont ete testes (11 6). les derniers 

codes de la fonction lineaire optimale D op constituent les codes de la fonction lineaire 

d'approximation (1 14). 

9. Methode selon la revendication 8 caracterisee en ce que la premiere phase comprend les etapes 
w modifiees suivantes : 

A1 - modifie I'etape A en operant une selection (102a) de deux couples de valeurs appartenant a 
ladite suite, tel qu'il existe au moins un couple additionnel intermediaire ayant une grandeur 
independante (Xi - Xe) comprise entre les grandeurs 
independantes dudit couple pour constituer au moins un triplet de couples de valeurs. 
'5 A2 - modifie I'etape A premierement en determinant (102b) une fonction lineaire annexe qui 

contient les deux couples selectionnes et deuxiemement en determinant des secondes 
erreurs entre les grandeurs dependantes des couples intermediates possibles et ladite 
fonction lineaire annexe : 

et Iorsque ces secondes erreurs sont toutes de meme signe, selection du couple 
20 intermediaire fournissant la plus grande seconde erreur, en valeur absolue, pour 

former un triplet de couples de valeurs forme du couple intermediaire et des deux 
couples selectionnes, 

. et Iorsque ces secondes erreurs sont de signes differents. reprise de la methode a 
I'etape A1, 

25 81 - I'etape B (104) est effectuee avec ledit triplet selectionne, 

C1 - modifie I'etape C (106) en selectionnant (106a) un couple additionnel dont la grandeur 
independante n'est pas comprise entre la grandeur independante des deux couples selec- 
tionnes. 

30 10. Methode selon la revendication 9 modifiee en ce que. prealablement a I'etape A1. lesdits couples de 
valeurs etant represents par des points (P1 - P6) dans un espace a deux dimensions, la premiere 
phase de la methode comprend une etape (100a) de determination d'une enveloppe inferieure et/ou 
d'une enveloppe superieure reunissant les points fes plus extremes de I'ensemble de points, la 
selection des deux couples de valeurs a I'etape A1 etant faite parmi les points adjacents appartenant a 

35 Tune ou I'autre enveloppe. 

11. Methode selon une des revendications 6 a 10 caracterisee en ce qu'a chaque valeur, prise dans la 
suite de valeurs de la grandeur independante, est associee un coefficient specifique de ponderation qui 
pondere lesdites premieres erreurs. 

40 

12. Methode selon une des revendications 6 a 11 caracterisee en ce qu'elle genere une fonction 
d'approximation comprenant une suite de fonctions lineaires de regression raccordees entre elles. 
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